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ABSTRACT 

We present spatially resolved observations of the canonical transition disk 
object TW Hya at 8.74 /im, 11.66 /im, and 18.30 /xm, obtained with the T-ReCS 
instrument on the Gemini telescope. These observations are a result of a novel 
observing mode at Gemini that enables speckle imaging. Using this technique, we 
image our target with short enough exposure times to achieve diffraction limited 
images. We use Fourier techniques to reduce our data, which allows high-precision 
calibration of the instrumental point spread function. Our observations span two 
epochs and we present evidence for temporal variability at 11.66 /zm in the disk 
of TW Hya. We show that previous models of TW Hya's disk from the literature 
are incompatible with our observations, and construct a model to explain the 
discrepancies. We detect marginal asymmetry in our data, most significantly at 
the shortest wavelengths. To explain our data, we require a model that includes 
an optically thin inner disk extending from 0.02 to 3.9 AU, an optically thick 
ring representing the outer disk wall at 3.9 AU and extending to 4.6 AU, and a 
hotter-than-disk-equilibrium source of emission located at ~3.5 AU. 
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Introduction 



Current planet formation theories are informed largely from observations of evolved 
planetary systems. Since the discovery of the first extrasolar planets (IMayor fc Quelozlll995l ) 
and the subsequent explosion of the field, theories of planet formation have grown and 
proliferated. Though the analysis of evolved systems is useful, perhaps the most effective 
way to test planetary formation theories is through the direct observation of young systems 
in the midst of planetary formation. 



Objects thought to be in this stage, so-called transitional disks (IStrom et al.l Il989l ). 
are intermediate between embedded pre-main-sequence stars and evolved planetary systems. 
These objects exhibit evidence of a low-density, optically thin gap in the innermost regions 
of their primordial circumstellar disks. The presence of this gap was originally inferred from 
the near- and mid- infrared (mid-IR) properties of the spectral energy distributions (SEDs). 
The SEDs of these objects exhibit a deficit of flux at wavelengths less than about 10 /im 
compared to the SEDs of classical T Tauri stars (CTTS), and a flux excess resembling CTT 
SEDs at wavelengths of ~20-60 fim. The deficit presumably occurs because the matter — and 
hence flux — in the inner regions of the disk is absent. The excess is caused by the outer disk, 
extending from the puffed up "edge" or "wall" at several AU to the extent of the disk at 
tens or hundreds of AU. 

Though their geometries were originally inferred from the properties of t heir unresolved 
SEDs , transition disk objects have since been resolved at various wavelengths. iThalmann et al 
(120 lOl ) analyze H- and i^-band data (1.65 /zm and 2.2 /zm, respectively) and resolve the gap 
in the transition disk LkCa 15. They confirm a disk with a ~46 AU t runcation radiu s , inside 
of which is a largely evacuated gap. This was also confirmed by Andrews et al.l (j2011al ) 
with imaging of LkCa 15 at 870 Atm. Several other transition disks have been resolved at 



2009, 2010; Andrews et al. 



2011b). 



millimeter wavelengths (Hughe s et al.ll2007t llsella et al.ll2008t iBrown et al.ll2009l ; llsella et al. 



The processes driving the formation of these gaps is not yet fully understood, though 
several hypotheses have been proposed (photoevaporation, grain growth, or a stellar or 
planetary companion). Many recent observations have been interpreted as evidence of a 
companion, supporting the hypothesis that young planets are important in the e volution 
and dissipation of the disk in these transition disk objects. iKraus fc Ireland! (120111 ) present 
the detection of a likely protoplanet located inside the known gap in LkCa 15. They use non- 
redundant aperture masking interferometry a t three epochs to disc over a faint companion 
located at ~ 16-21 AU from the primary star. iHuelamo et al.l (120111 ) detect a source in the 
L band (3.5 /xm) a t a separation of 6.7 AU — within the disk gap — from the transition disk 
T Chamaeleontis. lEisner et al.l (120091 ) present mid-IR observations of the transition disk 



-3 - 



SR 21. They spatially resolve the dust emission around this object, and suggest that the 
disk around SR 21 must be completely cleared within ~10 AU. They propose a disk with a 
large inner hole and a warm companion near the outer edge of the cleared region. 

TW Hydrae (TW Hya), a classical T Tauri Star, is one of the closest transition disks, 
located in the nearby (~50 p c: iMamaiek 120051) TW Hya Association. Due in part to its 
proximi ty and age C~1Q My r; I Webb et al.l 119991 ) . TW Hya has become an intensely studied 
object. ICalvet et al.l ( 120021 ) first proposed that TW Hya had a developing gap located at 
~4 astronomical units (AU), inferred from various features of the SED. More recently, 
spatially resolved observations of TW Hya have offered additional insight. The outer radius 
of the disk has b een resolved in mil l imete r wavelengths and determined to b e in the range 



of 70 - 140 AU flWilner et al.l I2000L 120031 : iQi et al.l 12004 : llsella et al.ll2009h . iHughes et al. 
(120071) resolved the inner edge of the gap in the circumstellar disk at 7 mm, corroborating 
Calvet et al.l (120021 ) 's determination of a hole at ~4 AU. Published c ontemporaneously w ith 



Hughes et al.l (120071 ) , iRatzka et al.l (120071 ) contradicted the claims of ICalvet et al.l ( 120021 ) in 
their analysis of interferometric observations of TW Hya in the mid-IR with the Very Large 
Telescope Interferometer ( VLTI) . They suggested that the disk gap occurs at a considerably 
smaller radius: 0.5 — 0.8 AU. At yet shorter wavelengths, Eisner et al.l (120061 ) measured the 
.fT-band (2 /im) visibility at a single baseline with the Keck Interferometer (KI), inferring an 
inner radius for the optically thin, evacuated region of 0.06 AU. 



Most recently, lAkeson et al.l (120111 ) presented near-infrared (near-IR) measurements 
from the Center for High Angular Resolution Astronomy (CHARA) array and the KI at 
various baselines. They model many of the past observations of TW Hya simultaneously 
with theirs. In combining the SEP and spatially resolved observations from the literature 
( lEisner et al.l 120061 : IHughes et al.l l2007t IRatzka et al.l 120071 ) , they review and extend past 
atte mpts at modeling t he circumstellar disk, producin g their own hybrid disk geometry. 
The lAkeson et al.l ( 120111 ) model is essentially the original ICalvet et al.l (120021 ) model with an 
added optically thick ring of emission at ~0.5 AU, at roughly the disk equilibrium temper- 
ature (see Section [3]). 

Modeling of spatially unresolved data for TW Hya (e.g. SEDs or higher resolution 
spectra) can yield very different conclusions. Spatially resolved observations are necessary 
to properly constrain the disk geometry. It is impor tant to note th e am biguity asso ciated 
with modeling source geometries from the SED alone: iBoss fc Yorkd (e.g.. Il993l Il996l ) found 
that the interpretation of infrared (IR) flux deficits as central gaps does not offer a unique 
solution, and that opacity and geometrical effects produce degenerate solutions fitting the 
SED of an unresolved system. We too find that this is true, and that the choice of dust 
species and thus inner disk opacity can have a significant effect on the determination of disk 
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geometry. Spatially resolved observations, ideally at multiple wavelengths, are needed to 
unambiguously constrain disk geometry. 

In order to further constrain models of TW Hya's circumstellar disk, and to resolve dis- 
crepant interpretations of observations, we present new speckle interferometric observations 
in the mid-IR. We observe TW Hya at three wavelengths (8.74 /mi, 11.66 /im, and 18.30 /im) 
and at two epochs (2007, 2009) using the Thermal-Region Camera Spectrograph (T-ReCS) 
instrument on the Gemini telescope. We resolve the disk at each wavelength, and construct a 
simple model to fit our observations. Combining our results with previous spatially resolved 
imaging and unresolved spectro-photometry, we constrain the properties of the inner disk 
(or lack thereof) in TW Hya. 

Drawing from t he literature, we firs t attempt to reprodu ce sim ple versions o f the m od- 
els first proposed bv ICalvet etaU hoO$ ): lUchida et all ( 120041 ) . and IRatzka et al.l (I2007I ). In 
the literature we identify two general classes of models: "Calvet-like" m odels with an opti 
cally thick wall at ~4 AU with a hole of optically thin material within (|Calvet et al.l 12002 



Uchida et al.ll2004j ; iHuehes et al.ll2007l : lAkeson et al-lbOlll iGorti et al.ll201lh : and "Ratzka- 
like" mo dels for which the r e exists a similar o pacity hole, but located much closer in at 
<1 AU (IRatzka et al.l 120071 ; lAkeson et al.l 1201 ll ). We approximately reproduce both model 
types, and compute synthetic mid-IR visibilities for comparison with our data; these are 
presented in Section 13.11 In Section 13. 2[ we describe a new model consistent with both the 
data from the literature and our new observations. 



2. Data 
2.1. Observations 

We present four observations of the transitional disk object TW Hya at three different 
wavelengths and two different epochs. A summary of these observations is presented in 
Table [TJ The wavelengths and epochs of the observations are: 8.74 /im (8.74) and 11.66 /im 
(11.66a), obtained in 2007; and 11.66 /im (11.66b) and 18.30 /im (18.30), obtained in 2009. 
For each observation, we also took off-target nod pointings for infrared sky background 
subtraction. Each nod, or pointing, consists of several short exposures ("frames") that are 
combined to yield a high resolution speckle image. It is the individual frames, statistically 
combined, that are ultimately used in our analysis. For each dataset, we observed a bright, 
point source calibrator in order to calibrate and remove instrumental and atmospheric effects. 

We collected our data using a custom observing mode on the T-ReCS instrument at 
Gemini. We used short integration times (tj nt ~172 msec) to freeze the motion of the 
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atmosphere and thus achieve high spatial resolution, enabled by the diffraction limit of the 
telescope. We employed a dither / nod pattern that moved observed objects around four 
positions on the detector. For the 18.30 /xm dataset, we employed standard chopping, and 
our frame integration times were a factor of ~10 longer than at the other wavelengths. Total 
integration times are included in Table [TJ 

The fluxes for our SEP comparisons were obtained from the Spitzer Space Telescope 
Infrared Spectrograph, from lUchida et al.l (120041 ) (see their Figure 2). For our calculations 
of our best-fit model, we increase their stated errors of ~10% to ~30% to account for the 
larger number of points in the SED as compared to our resolved dataset, and to weight higher 
the value of the resolved data oy er the unresolved data in our fits. The errors displayed in the 
Figures are lUchida et al.l (120041 ) 's ~10%. We also use mid-IR interferometric measurements 
for our anal ysis. These da t a wer e obtained with the MIDI instrument at the VLT, and were 
presented in iRatzka et al.l (120071 ) . 



2.2. Reduction and Analysis 

2.2.1. Fourier Analysis 

We reduce our data using Fourier analysis techniques. Each individual frame has an 
exposure time that alone is too short to provide a significant signal to noise ratio. However, 
a single long exposure would yield an image with insufficient angular resolution due to the 
effects of the Earth's turbulent atmosphere. One cannot naively add together the short, 
individual frames, as the target centroid shifts on the sky due to atmospheric turbulence. 
We instead combine the power spectra of individual frames. This method of addition is 
independent of translations of the image centroid position. Furthermore, analyzing the data 
in Fourier space allows us to remove the instrumental point spread function (PSF) with a 
point-source calibrator with high accuracy. In Fourier space, we divide the power spectrum 
of the source by the power spectrum of a point source calibrator. This is considerably simpler 
than deconvolution in the image plane. 

The first step in our data reduction procedure is to remove the sky background. We 
obtained our data in such a way that each series of frames (short exposure images) of a 
target is paired with a slightly offset series of frames of the same target; we call these two 
sets "adjacent nods". This offset causes the target to appear on two different regions of 
our detector. We use these adjacent nod pairs to overcome the high sky background in 
the mid-IR by subtracting from each target its adjacent nod. As each nod contains many 
frames, the median value of all the frames in a particular nod is used for the background 
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subtraction. The median value of an adjacent nod is subtracted from each individual short- 
exposure frame. For example, in a single nod pair, the median value for all the frames in 
Nod Type "A" is subtracted from each frame in Nod Type "B", and likewise the median 
value of "B" frames is subtracted from frames in "A" . This process is repeated for each pair 
of nods. 

Next, we flag unusable frames, without actually modifying any of the data used for 
analysis. We first attempt to locate a bright point source in the median image of a pointing: 
either the target — TW Hya — or a calibrator. To locate the position of the star in a pointing 
(or nod), we calculate the median image for that nod, subtract from the median image the 
mean value of all the pixels, and set negative pixel values to zero. We identify "hot pixels" 
as pixels with values more than three standard deviations larger than the mean pixel value 
of the resultant image. These "hot pixels" are assigned a value equal to the average value 
of their nearest neighbors using image convolution. This technique removes isolated "hot 
pixels" while leaving signal from a star relatively unaffected. The position of the star is then 
set to the location of the maximally valued pixel in the image. We classify the pointing as 
unusable if, after this process, there exist no pixels that are >5-cr deviations from the mean, 
or if the location of the point source determined by our algorithm is closer than the width of 
our subimages to the detector edge (see next paragraph). If we fail to detect a bright point 
source in a particular pointing, we flag all the frames in that pointing as unusable. This 
could happen due to poor seeing in these exposures, the presence of clouds, or some other 
effect. We note again that the steps described in this paragraph do not ultimately affect the 
data used for analysis, but only for data flagging. After a star is located, we return to the 
background-subtracted data described at the end of the previous paragraph. 

After eliminating flagged data, we are left with frames containing usable point sources. 
For each frame, we cut out a "postage-stamp" size subimage from the raw telescope image, 
centered on the point source. These subimages are 64x64 pixels, or 5.76x5.76 arcseconds at 
T-ReCS's plate scale, or ~309x309 AU at the distance of TW Hya. We are then left with 
order hundreds (ranging from 24 to 3400) of short-exposure frames for target and calibrator 
for each dataset, cropped and centered at the star's position. 

For each frame, we subtract a residual sky value in addition to the initial IR sky back- 
ground subtraction: the median of pixel values in the outer regions (48 pixels or farther 
from the image center of each subimage). We then apply a Hanning window to the image 
by multiplying by the two-dimensional Hanning function: 



where D is the size (diameter) and (x ,y ) is the center of the Hanning window. For our 




(1) 
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analysis, D — 16 pixels (roughly 1.4 arcseconds or 77 AU). The Hanning window has some 
useful properties: it is unity at the center and falls smoothly to zero at the edges. This 
step removes spurious Fourier signals by smoothing sharp brightness transitions. We take 
the Fourier transform of the Hanning windowed image, and record the squared amplitude at 
each pixel (two dimensional power spectrum, or power image). We perform steps identical 
to those described above for an adjacent, equally sized, background-subtracted subimage of 
blank sky from the same parent image to produce a sky power image. 

After computing the Fourier amplitudes, we use sigma-clipping to discard discrepant 
data. For each dataset, we calculate the mean and standard deviation of the squared visibility 
as a function of baseline for both target and calibrator. If any individual frame deviates >3-a 
from the mean at any moderate, well-behaved baseline (between m and 5 m), we exclude 
it from the sample. This step excludes a small fraction of our resulting frames: usually zero, 
but occasionally as high as 2%. 



2.2.2. Calibration 



After the reduction procedure outlined above, we are left with several hundred power 
images, for both target and calibrator. Additionally, we have power images for adjacent, 
blank sections of sky from parent images for each of these postage-stamp images. We first 
remove any bias due to detector artifacts by subtracting the sky power image from the 
object power image. We then subtract a residual bias: any non-zero value for the visibility 
amplitude at baselines well beyond our sensitivity (9 m and longer). We set negative values 
of the target's power images to zero, and negative values of the calibrator to a small value 
(10~ 5 ; to avoid eventual division-by-zero computational errors). We then normalize the 
target and calibrator power images by the maximum value in each respective set. 

We calculate an azimuthal average of each power image by computing the mean of the 
pixels in an annulus corresponding to a particular baseline length. The power obtained in 
each annulus becomes the visibility amplitude at a corresponding baseline. We do this to 
increase the signal to noise and to make our plots more straight-forward to interpret. Further- 
more, to rough approximation, it is an acceptable assumption that our target — a star with 
a near ly face-on circumstellar disk — is indeed symmetric about the azimuth. ICalvet et al. 
( 120021 ) also assume a face-on orientation, in agreement with estimates of TW H ya's inclina- 
tion in the literature ( IKrist et al.ll2000l : iQi et al.ll2004t iPontoppidan et al.ll2008l ). We discuss 
possible departures from our assumption of azimuthal symmetry in Section 13.31 



We use the azimuthally averaged visibility curves (one for each frame) to calculate 
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a mean and standard deviation of the mean (SDOM, or standard error) for target and 
calibrator. We then use these derived values (V trg , cr trg , V ca \, <r ca i) to calculate the calibrated 
visiblity (Vt rg /Kai) an d its associated statistical uncertainty. The SDOM is calculated by 
dividing the usual standard deviation by the square root of the number of frames (vjV^). 

The total statistical uncertainty associated with the calibrated visibility is given by the 
propagation of the two SDOMs. We calculate the propagated error using the usual relation, 



for f = f(a, b) and for uncorrelated errors in a and b. 

We also estimate a systematic error for our experimental setup, by comparing calibrator 
observations. For example, variation in seeing or sky background structure over timescales 
longer than our source-calibrator cadence will lead to systematic errors. Similar to the 
procedure described above, we average the azimuthal average of power images according 
to telescope pointing (i.e., we group calibrator frames of a single pointing together). We 
divide pointing-averaged calibrator visibility curves by the mean visibility curve of all the 
calibrators. We use the SDOM — not the standard deviation — of the several calibrator point- 
ings as a systematic error: the magnitude at which we expect an uncertainty to exist for a 
given set of observations. For each baseline, the visibilities of the calibrators do appear to 
be normally distributed. We add this error in quadrature with the statistical uncertainties 
described above. The final error in our measurements is dominated by these estimated sys- 
tematic errors. We present the data used for this method of systematic error estimation in 
Figure [H This plot shows the visibility amplitude of each calibrator pointing, divided by the 
mean visibility amplitude of all the calibrators, as a function of baseline, for each dataset. 

We expect the calibrator errors to appear correlated to some extent. That is, if a 
calibrator has a larger visibility amplitude at one baseline, it likely has a larger visibility at 
other baselines as well, as seeing variations are an important cause of the dispersion. One 
other possible cause for correlated errors is our application of a Hanning windown in the 
data reduction process. To check this effect, we perform our analysis without the Hanning 
window; we notice a slight decrease in apparent correlation. The size of the calculated error 
bars did not change appreciably, however, and we keep the Hanning window in our reduction 
process. 




(2) 
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Fig. 1. — In this Figure, we illustrate the method used to estimate systematic errors in our 
visibility measurements. For each observation block, we plot the value of the visibility for 
a calibrator as a function of baseline, divided by the mean value for all the blocks. We 
perform this test for each dataset, shown in the separate panels: 8.74 (blue, top left), 11.66a 
(black, top right), 11.66b (black, bottom right), and 18.30 (red, bottom left). The value of the 
systematic error estimate is the standard error (or standard deviation of the mean, SDOM) 
of the observation blocks at each baseline. 
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2.3. Size of Emitting Region in TW Hya 

We estimate the size of the emitting region at each wavelength by fitting a ring of fixed 
thickness (width to radius ratio, / =0.11) to each visibility dataset (see Equation [6]). We 
performed a chi squared minimization where the size of the ring of emission was the only free 
parameter. Error bars (l-cr) were derived by finding the ring size value that corresponded 
to the minimum chi squared plus one. The results of our fit can be seen in Table [2] and a 
comparison of our data with the best fit ring visiblity curves can be found in Figure |2j 



2.3.1. Variability 

Since we have two epochs at 11.66 /im, we are able to constrain variability over a 
timescale of roughly two years (~720 days). Variability is observed between these two epochs 
at the >2.5-cr level (Figure EJ Table [2]). For simplicity, we use both epochs of 11.66 //m data 
in our modeling. We speculate as to the cause of this variability in Section HI 



3. Modeling 



We model the TW Hya system using largely geometric models. We approximate the 
central star as the Kurucz stellar atmosphere of a K7 star with radius = R Q , temperature 
T* = 4000 K, surface gr avity log g (cm s -2 ) = 4.5, at a distance d = 53.7 pc (jWebb et al. 



1999 



van Leeuwenl 120071 ). We note, as has been mentioned previously in the literature (e.g., 



Sitko et al.l 120000 . that TW Hya is a known variable star, so we do not necessarily expect 
a Kurucz model to exactly model the emission of the star. Also, the shorter-wave flux 
density measurements from the literature were not taken contemporaneously with the mid- 
IR spectrum that comprises the bulk of the disk emission. While this Kurucz model does 
not perfectly reproduce the stellar flux at all wavelengths, and different analyses of TW Hya 
use slightly different stellar model parameters, it is sufficient for our modeling purposes. 

Our models of TW Hya are composed of an optically thin disk (approximated as a series 
of concentric rings which follow a known temperature-radius relationship) and an optically 
thick ring of a single temperature representing the directly illuminated edge of the optically 
thick disk. For the optically thin component, we must choose and apply dust opacities. We 
discuss the details of this choice below. 
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Table 1. 


Observations 


of TW Hya 






Date 


Calibrator 


A (/im) 


tint (msec) 


N nod 


N bm 




71 nt (sec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


9 May 2007 


II Hya 


8.74 


173 


18 


702 


351 


61 


9 May 2007 


II Hya 


11.66 


173 


30 


1170 


1131 


195 


9,18 April 2009 


HD 92036 


11.66 


181 


100 


3400 


2505 


454 


22,23 May 2009 


HD 92036 


18.30 


1813 


76 


228 


228 


413 



Note. — The data used in our analysis. We observed the target, TW Hya, at three different 
wavelengths and at two different epochs. We discard unusable frames, according to criteria in 
Section 12.21 Columns are: (1) date of observation; (2) name of calibrator star; (3) wavelength 
of observation in microns; (4) integration time of an individual frame in milliseconds (msec); 
(5) number of nods in the observation of the target; (6) number of individual frames in the 
observation of the target; (7) number of frames used for analysis, after removing flagged frames; 
and (8) total integration time for target using all usable frames, in seconds (sec). 



Table 2. Single Wavelength Ring Fit Results 



Observation Name 


Rin [AU] 


/ 


(1) 


(2) 


(3) 


871 


1.98 ±0.04 


0.11 


1166a 


2.56 ±0.23 


0.11 


1166b 


0.71 ± 0.50 


0.11 


1830 


3.35 ± 0.31 


0.11 



Note. — The results of a single wave- 
length ring fit to our observations. Error 
bars are derived by finding the value (in ra- 
dius) on the chi squared contour that corre- 
sponds to one plus the minimum chi squared 
value. Columns are: (1) name of observa- 
tion; (2) best fit radius with 1-<t error bars; 
(3) the ratio of ring width to inner radius, 
set to a constant for these illustrative pur- 
poses. We note the >2.5-<r difference in size 
between our two epochs of 11.66 fim obser- 
vations. 
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Fig. 2. — This figure shows the visibilities predicted by ring models of various sizes. The 
rings emit at a single wavelength and have width to radius ratio of / = 0.11. This exercise 
shows that our 8.74 /xm observations (blue, top left) are significantly more resolved than all 
of our other epochs The four separate observations are: 8.74 (blue, top left), 11.66a (black, 
top right), 11.66b (black, bottom right), and 18.30 (red, bottom left). 
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3.1. Existing Models 



We first attempt to reproduce the large cavity model of lCalvet et al.l ( 120021 ) ; lUchida et al. 



( 12004 ) . We use geometric mode l parameters and du st opacities from their work. We extract 



dust opt ical depth dir e ctly fr om lUchida et al. d200# l. which presents results consistent with 
those of ICalvet et al.l (120021 ). This model consists of emission from three components: 1) 
a star, 2) an optically thin disk, extending from the dust sublimation radius at 0.02 AU 
to the directly illuminated disk wall at 3.3 AU, and 3) an optically thick rin g at 3.3 AU, 



(2002 


); 


Uchida et al. 


(2004) 



by the Kurucz stellar model to vary so that it matches the observed flux den sity at 4 — 5 /zm; 



this is similar to their own strategy and the procedure of ISitko et al.l (120001 ) 



The flux from the optically thick disk wall is approximated by an unmodified blackbody. 
The flux density of an annulus of thickness dR is given by: 



2ir 

dF ^ = rx^B x (T dust )RdR. 



(3) 



In the case of the optically thick ring, T\ — l for all values of A. For the optically thin disk, 
the temperature profile is given by 



. R* 

1 disk — J- * 1 7^ 



1/2 



(4) 



where T* is the temperature of the star, i?* is the radius of the star, and R is stellocentric 
radius. To calculate the flux from the disk, we use small concentric annuli. The temperature 
structure at every point in the optically thin disk behaves according to Equation HI and the 
fluxes from each annulus are summed together to obtain the total disk flux. The temperature 
for the directly illuminated rim or edge of the optically thick disk is given by 



T ■ 



1 R 



1/2 



(5) 



We ar e motivated in the cho ic e of these temperatu re profiles by the methods of lChiang fc Goldreich 
( 119971 ): ICalvet et al.l ( 120021 ): lEisner et al.l (120061 ) (see Equation 1, Equations 1 and 3, and 
Equation 1, respectively). 



1 To obtain the optical depths used bv lUchida et al.l (|2004l) . we first record the reported flux density from 
their optically thin disk component (see their Figur e 2, top panel ) . We then use these published flux densities, 
an assumption of constant surface density (as in ICalvet et al.l rt2002h l. and their stated disk geometry, to 
extract values for the dimensionless wavelength-dependent disk optical depth, t\, that they use. 
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All model components can be viewed as collections of single-temperature rings that obey 
the model temperature profile. Model fluxes are computed by summi ng Equation p ove r all 
annuli. The squared visibility for a ring (or annulus) is given by (see Eisner et al.ll2003l ): 

2 



V T 



ring 



x [(1 + /) Ji [(1 + f7T6 in r uv ] - Ji(7rt? in r ra )] 



(6) 



7rr WJ in (2/ + /2) 

where r uv = yju 2 + v 2 = B ■ s/X is the u-v radius, A is the wavelength, 9 is the angular 
diameter of the object in radians, / defines the ratio of radial thickness to size of the annulus 
(/ — W/R, where W is the width of the annulus and R is the inner radius of the annulus), 
J i is the Bessel function of the first kind of order one. The normalized visibility for the 
entire model is the flux-weighted average of the visibilities produced by each annulus. 
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( 12004] ) ). we generate a synthetic SED and visibilities (Figure [3]). The large bottom panel 
shows the SED, indicating the stellar flux density {thin solid line), the optically thin, low 
density inner disk (dotted- dashed line), the outer disk wall (dotted line), and the total flux 
(thick solid line). The flux data are also shown (blue circles). In the top half of the figure, we 
show our four separate observations: 8.74 (blue, top left), 11.66a (black, top right), 11.66b 
(black, bottom right), and 18.30 (red, bottom left). In each, the flux weighted visibilities of 
the different model components are indicated by different line types, as in the SED plot. The 
insets in the top right of each panel show zoom-ins of the data presented in this work, while 
the long baseline points at ~45 m are from lRatzka et all (120071 ). Vertical, dashed lines in the 
SED panel show the wavelengths of our observations at 8.74 /mi, 11.66 /im, and 18.30 /im. 



Our estimate of ICalvet et all (120021 ) ? s model reproduces the SED as well as done by the 
original authors. This is not unexpected, as we obtained the opacities, physical geometries, 
and temperatures directly from their work. More interesting are the mid-IR visibilities 
produced by this model, shown in the top half of the figure. This model reproduces well 
the 18.30 /xm data, a tracer of cooler portions of the disk, dominated by emission of the 
transition disk wall. Similarly, this model reproduces rather well both epochs of our 11.66 /im 
data. One epoch (1166b) is reproduced better than the other (1166a), but we note that 
these two datasets are in fact inconsistent with each other due to temporal variability (see 
Section r2.3.1l) . This model is clearly inconsistent, however, with our very resolved 8.74 /im 
observations, as seen in the top left panel of Figure [31 



3.1.1. Small Cavity Models 



Though subsequent work has c onfirm ed the validity of ICalvet et all (120021 ) 's model (e.g., 
Hughes et al.ll2007l ). iRatzka et al.1 (120071 ) claimed that the evacuated cavity was much closer 
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Fig. 3.— The calculated SED and visibility curves for our realization of ICalvet et al. 



(j2002l )'s model. See Table [3] for model parameters. The large bottom panel shows the SED, 
indicating the flux from the star (thin solid line), the optically thin inner disk (dotted line), 
the optically thick, directly illuminated disk wall (dashed line), and the total (thick solid 
line). The flux density data are also shown (blue circles). In the top half of the figure, we 
show our four separate observations: 8.74 (blue, top left), 11.66a (black, top right), 11.66b 
(black, bottom right), and 18.30 (red, bottom left). In each, the flux weighted visibilities of 
the different model components are indicated by different line types, as in the SED plot. 
The insets in the top right of each panel show the data presented in this work, while the 
long baseline point at ~45 m is from lRatzka et al.l (120071 ). Vertical, dashed lines in the SED 
panel show our wavelengths of observation. 
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to the star, at <1 AU. Since iRatzka et al.l ( 120071 ) 's model was based largely on mid-IR 
flux measure ments, and presented resolved mid-IR data, we compare their work against our 
observations. IRatzka et al.l (120071 ) used a more complicated disk modeling technique than did 



Calvet et al.l ( 120021 ); we use our simple geometric models to distill out the most important 



elements. The main sources of the mid-IR emission in the model of IRatzka et al.l (120071 ) 
were the optically thick, directly illuminated disk wall, and its optically thin atmosphere. 
We were thus motivated to produce emission at <1 AU from nearly coincident components, 
one optically thick and one optically thin. 



Our methods to reproduce the IRatzka et al.l (120071 ) model are similar to the process 
described previous ly in Section 13.11 w ith one exception. As the bulk of the optically thin 
emission from the IRatzka et al.l (120071 ) model originates from the disk atmosphere, we use 



an optically thin ring at a single temperature (-R a tm ; ^atm) instead of the optically thin inner 
disk described above (Ri n , T in ). The temperature of this ring is a model parameter and is 
free to vary. We also vary the surface density of the optically thin ring, E atm . We vary 
these free parameters and minimize chi squared of the model c ompared to the mid-I R flux of 
TW Hya and the long-baseline (~45 m) v isibility points from IRatzka et al.l (120071 ) . As our 
motivation was to reproduce the model of IRatzka et al.l (120071 ). we do not include our new 
visibility data in the minimization of chi squared. 



Ratzka et al.l (120071 ) used different optical depths than did lCalvet et al.l ( 120021 ) ; IUchida et al 
( 120041 ). and we obtained and used these in our modeling (T. Ratzka 2011, private commu- 
nication). As in the modeling described above, we allow the flux density produced by the 
Kurucz stellar model to vary so that it matches the observed flux dens ity at 4 — 5 /im; this 
is similar to their own strategy and the procedure of ISitko et al.l (120001 ) . 



We performed a chi squared minimizatio n and allowed the fo llowing; parameters to vary 
freely to generate a model consistent with the IRatzka et al.l (120071 ) geometry: Kurucz scaling, 
-Ratm, T atm , S atm , -Rwaii, /wail, ^waii- These properties and their best-fit values are listed in 
Table EJ 

We present the calculated SED and visibility curves for our realization of this model in 
Figure HI The large bottom panel shows the SED, indicating the stellar flux density (thin 
solid line), the optically thin, transitional disk wall atmosphere (dotted- dashed line), the 
transitional disk wall (dotted line), and the total flux (thick solid line). The flux density 
data are also shown (blue circles). In the top half of the figure, we show our four separate 
observations: 8.74 (blue, top left), 11.66a (black, top right), 11.66b (black, bottom right), 
and 18.30 (red, bottom left). In each, the flux weighted visibilities of the different model 
components are indicated by different line types, as in the SED plot. The insets in the top 
right of each panel show the data presented in this work, while the long baseline points at 
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~45 m are obtained from lRatzka et all (120071 ). Vertical, dashed lines in the SED panel show 
our wavelengths of observation. 



Our realization of IRatzka et all ( 120071 ) 's does not perfectly match the model SED pre- 
sented in iRatzka et all ( 120071 ) (particularly at 16 — 18 jttm), but does a reasonably good job 
considering the complex disk properties (e.g., dus t settling) we did no t include. Our model 
visibilities do reproduce well those presented in IRatzka et all ( 120071 ) at long baselines at 
8.74 /im and 1 1.66 u rn (see their Figure 8 for comparison). As can be seen in this figure, the 
Ratzka et all (120071 ) data show t he emission at ~8.7 4 /im to be more resolved than that at 
~1 1.66 /mi, though the model of IRatzka et all (120071 ) does not reproduce this behavior. We 
see this trend of very- resolved 8.74 /mi emission in our own observations as well. 

It is clear, however, that this small cavity model is inconsistent with nearly all our data. 
Though it is consistent with one epoch of our 11.66 /im data (1166b), it is very inconsis- 
tent with all the other observations. It is unresolved at all baselines and all wavelengths 
relevant to this work. While we did not include an outer disk, this does not impact the 
observed discrepancy because the emission from these cool, outer regions is not relevant for 
the wavelength ranges we are interested in. This omission may have slightly worsened our 
reproduction of the small cavity model's 18.30 /im emission, but outer disk components at 
cooler temperatures will not cause the model to show an 8.74 /im component as resolved as 
our data show. 

We suspect that the difference in dust o pacity choice s is th e principal reason that two 
distinct disk geo metries — the ^4 AU holes in lCalvet et all ( 120021 ) models, versus the <1 AU 



hole inferred by IRatzka et al.l (120071) — can reproduc e the same SED. I n Figure El we com 



pare the optical dep ths used by IRatzka et al.l ( 120071 ) to those used by ICalvet et al.l ( 120021 ) ; 
Uchida et al.l (120041 ) . In order to examine the difference s between the opacit y choices in these 
two works, we have converted the optical depths from IRatzka et al.l (120071 ) to a dimension- 
less optical depth with an assumption of surface density (informe d from their work) . We 
then scaled t h is val ue so that the mean values of optical depth for ICalvet et al.l (|2002l) and 
Ratzka et al.l ( 120071 ) match, thus facilitating comparison. The IRatzka et al.l (120071 ) optical 
depths show a relative paucity around 11 — 12 /mi that explains why a smaller (<1 AU versus 
~4 AU for a large cavi ty model) transition a l disk rim is used in that model. We note that the 
optical depths used by lCalvet et al.l (120021 ) ; IUchida et al.l (120041) were der ived simultaneously 
with the other disk parameters, while those used by IRatzka et aL ( 2007 ) were not. 
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Fig. 4— The calculated SED and visibility curves for our realization of iRatzka et al. 



(120071 ) 's model. See Table [3] for model parameters. The large bottom panel shows the SED, 
indicating the flux from the star (thin solid line), the optically thin disk wall atmosphere 
(dotted line), the optically thick, directly illuminated wall (dashed line), and the total (thick 
solid line). The flux density data are also shown (blue circles). In the top half of the figure, 
we show our four separate observations: 8.74 (blue, top left), 11.66a (black, top right), 11.66b 
(black, bottom right), and 18.30 (red, bottom left). In each, the flux weighted visibilities of 
the different model components are indicated by different line types, as in the SED plot. 
The insets in the top right of each panel show the data presented in this work, while the 
long baseline point at ~45 m is from lRatzka et all (120071 ) . Vertical, dashed lines in the SED 
panel show our wavelengths of observation. 
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Fig. 5. — The opt i cal depth used by t he works that initially pr esented large centra l cavity 
(ICalvet et al.ll2002[ lUchida et al.ll2004j ) and small central cavity ( iRatzka et al.l 120071 ) models 
to explain the mid-IR emission from the TW Hya disk system. We note the discrepancies 
in choice of optical depth by these authors. As the optical depth and disk geometry are 
degenerate in model-fitting an unresolved SED, these two geometrically distinct models are 
both able to reproduce the SED. To examine t h e rela tive differences between choices of 
optical depth, we have scaled the iRatzka et al.l ( 120071 ) optical depths so that the mean 
values of each optical depth plot match. 
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3.2. New Disk Model 



Neither of the previous models considered (ICalvet et al.ll2002t IRatzka et al.l 120071 ) ad- 
equately fit our new observations. The hybrid model of Akeson et al.l (2011) also cann ot 
explain our very resolved 8.74 /im data. This is not surprising, as lAkeson et al.l (120111) 's 



Calvet et al. 



(2002 


); 


Ratzka et al. 


(2007) 



emission at 8.74 fim, without producing heavily resolved emission at larger wavelengths, we 
need a hot component at large stellocentric radius. We will show that this "hot" component 
must be considerably hotter than predicted by equilibrium disk models. 



We adopt the optical depths from lUchida et al.l (120041 ) (see Section [37Tj) . We assume 
a Kurucz stellar atmosphere with the scaling left as a free parameter. In addition to the 
components described in Section 13. 1[ we include an additional optically thick ring (-Rhot and 
^hot)- We also vary the surface de nsity of the optically thin inner disk, £;„ , but leave the 



dust sublimation radius R in used by lCalvet et al.l ( 120021 ); lUchida et al.l (120041 ) unchanged, at 
0.02 AU. We allowed the following parameters to vary freely (between physically reasonable 
limits) to generate our best-fit model: Kurucz Scaling, R wa xi, /wall, T wa a, S in , -Rhot, /hot, and 
Thot- For definitions of these parameters, see Table [3j 

After a thorough exploration of parameter space, we found a model that explains the 
existing data. The observables associated with this model are presented in Figure [6j In this 
figure we show the calculated SED and visibilit y curves for the mod el presented in this work. 



This mo del is a large cavity model, similar to ICalvet et al.l ( 120021 ) and consistent with the 



results of iHughes et al.l (120071 ). The directly il luminated disk wall is represented by a ring of 



optically thick emission at 3.9 AU, and, as in lCalvet et al.l (120021 ). there is an optically thin 



disk of emission at low surface densities that extends from the dust sublimation radius to the 
outer disk wall (0.02 — 3.9 AU). The crucial addition to this model is an additional optically 
thick ring of emission inside the disk wall, at ~2.9 AU. In order to fit our 8.74 fim data, 
this ring is at a hotter than equilibrium temperature (which would be ~150 K, compared 
to ~570 K for the hot ring). We also find that the best fit is achieved with an unsealed 
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); 
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(2007) 



We note that with this model we produce a good fit to the SED, discrepant only in the 
6-8 /j,m region where possible stellar variability is a larger co ncern. Even so , the d iscrepancy 
in the SED fit is comparable in magnitude to that of the lUchida et al.l ( 120041 ) deviation 
from observations in the ~13 jim region. We not e that our fit to the long baseline mid- 
IR visibility amplitudes fr o m IRatzka et al.l (120071 ) is comparable to the model visibilities 
produced by IRatzka et al.l (j2007l )'s own model. Most notably, however, we point out the 
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simultaneous high quality fit of the long baseline mid-IR visibility amplitudes with the short 
baseline mid-IR visibility amplitudes presented in this work. Again, this was accomplished 
by forcing the ring to be hot enough to peak at shorter wavelengths — closer to 8.74 fim — 
allowing the short bas eline 8.74 /jm emiss ion to be resolved, as we see in both our data 
and the VLTI data of iRatzka et all (120071 ). Additionally, despite the increased degrees of 
freedom, the reduced chi squared value for our new disk model is ~1% better than the value 
we calculate for the Calvet model. The reduced chi squared value for the Ratzka model is 
considerably larger. The value of the chi squared is largely driven by the unresolved epoch of 
our 11.66 /iin micron data, however. The less-resolved Calvet model fits this epoch slightly 
better. If we fit the average of the two 11 micron epochs, instead of attempting to fit both 
simultaneously, our new model has a reduced chi squared that is a factor of 3 lower. 



3.3. Disk + Companion Model 

The models described thus far are all symmetric about the azimuth. As described in 
Section 12.2.21 the reason we made this assumption was to decrease the uncertainty associ- 
ated with our measurements of baseline-dependent visibility and to make our results more 
straight-forward to interpret. However, our solution of a hotter-than-equilibrium, geometri- 
cally thin ring of emission at a radius just inside the directly illuminated disk wall does not 
seem physical. That is, we cannot think of a physical process that could heat an annulus of 
matter far beyond the temperature achieved for a disk in equilibrium with stellar emission. 

A self-luminous companion, on the other hand, can lead to such heating. Such an 
object would clearly not be azimuthally symmetric. The presence of a companion should 
thus manifest itself in our data through departures from azimuthal symmetry as a function 
of wavelength. 

If all of our model components except the companion are azimuthally symmetric, we 
can predict relative degrees of asymmetry in the different wavelengths of observation. We 
expect the azimuthal asymmetry to be largest at 8.74 /im; at 11.66 and 18.30 /im, the flux 
contribution from the hot inner component is small compared to fluxes from other, symmetric 
model components (see Figure [6]). 

To test this, we create a model with a companion that has the same temperature and 
surface area as our added ring component and thus the disk plus companion model produces 
an SED identical to the model described in Section 13.21 For this reason we can utilize only 
our visibility data to constrain properties of the companion's stellocentric radius and position 
angle (PA); the SED will not change as these parameters are varied. 
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Fig. 6. — The calculated SED and visibility curves for the new model presented in this work. 
See Table [3] for model parameters. The large bottom panel shows the SED, indicating the 
flux from the star (thin solid line), the optically thin inner disk (dotted line), the optically 
thick, directly illuminated wall (dashed line), the hot, optically thick ring (dashed- dotted 
line), and the total (thick solid line). The flux density data are also shown (blue circles). 
In the top half of the figure, we show our four separate observations: 8.74 (blue, top left), 
11.66a (black, top right), 11.66b (black, bottom right), and 18.30 (red, bottom left). In each, 
the flux weighted visibilities of the different model components are indicated by different line 
types, as in the SED plot. The insets in the top right of each panel show the data presented 
in this work, while the long baseline point at ~45 m is from iRatzka et al.l ( 120071 ). Vertical, 
dashed lines in the SED panel show our wavelengths of observation. 
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We then calculate the degree of asymmetry generated by the presence of a companion 
in this model. Similarly, we measure the asymmetry in our reduced data. We follow the 
procedures outlined in Section 12.21 with one exception. Instead of performing an azimuthal 
average of the power image, we calculate the radial average as before, but include only small 
regions (of width 45 degrees) of azimuth separated by 23 degrees. That is, we calculate 
the radial average (power as a function of baseline) in different "wedges" of azimuth. We 
perform these steps for both a synthetic image of the disk plus companion model, and for 
our data. 

Using the azimuth wedges for both synthetic model and reduced data, we compare the 
disk plus companion model to the data. We fix all the components of the model except the 
companion stellocentric radius and position angle. We then perform a best-fit chi squared 
minimzation over these two freely varying parameters (PA and stellocentric radius) for all 
epochs of our data (see Tabled]). We find a best-fit companion radius of ~3.5 AU, and a 
PA of ~90 degrees. We note that as we have limited phase information using this technique, 
that we are unable to distinguish between PAs separated by 180 degrees. Thus, our best-fit 
PA is ~90 degrees or ~270 degrees. Position angle is measured counter-clockwise eastwards 
from north. 

Since our data span two epochs separated by ~720 days, we do not fit either epoch of the 
11.66 /j,m data perfectly. Further, if a companion indeed explains our data, then this object 
would undergo orbital evolution between our two epochs. For these reasons, we additionally 
fit PA separately for each epoch for a fixed radius: 3.5 AU, the best-fit companion radius 
for all the data. This exercise yields a best-fit PA of 90 degrees for our 2007 data (874 and 
1166a), but the PA is not strongly constrained in the 2009 data (1166b and 1830). This is 
unsurprising, as the 2009 epoch does not contain short-wavelength observations in which the 
flux from the hot companion would be significant. 

To more clearly illustrate asymmetry in our multi-epoch data and to ameliorate ap- 
parent discrepancies due to fitting our disk model to the combined 2007 and 2009 datasets, 
we perform a simple "de-trending" of the model visibility amplitudes. That is, we force the 
synthesized visibilities to match the shape of the visibilities obtained through our data reduc- 
tion by applying a multiplicative scaling at each baseline. As visibilities greater than one are 
unphysical, we set detrended model values that exceed one to be equal to one. In Figure [7J 
we show the expected magnitude of the asymmetry and compare this to the magnitude of 
the asymmetry in our observations. 

In this figure, the points with error bars show the degree of asymmetry revealed in 
our data reduction, using the wedge comparison technique described above. We choose the 
best-fit PA for each epoch for one wedge, and the best-fit PA plus 90 degrees for the second 
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Fig. 7. — This figure shows the asymmetry in our observations. The two sets of points with 
error bars in the figure show the normalized, calibrated visibility amplitude as a function 
of baseline, as in Figures [3j HJ and |6j Here, though, instead of azimuthally averaging as 
described in Section I2.2[ we perform a radial average along two narrow slices of angle in 
azimuth, offset by 90 degrees. Each curve thus shows how resolved the image is along each 
direction in the Fourier power image. This figure shows that the 8.74 fim emission is very 
resolved along one direction, and quite unresolved along another orthogonal direction. These 
data are inconsistent with an azimuthally symmetric ring of emission. The solid and dotted 
lines show the "de-trended" (see Section 13. 3j) predicted degree of asymmetry from a best-fit 
disk plus companion model. 
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wedge. Though not highly statistically significant, this figure shows that the asymmetries are 
consistent with our expectation from an asymmetric source of emission at short wavelengths: 
that the asymmetry in the 8.74 /xm data is greatest, smaller at 11.66 /xm, and yet smaller at 
18.30 fim. The solid and dotted lines are the visibilities obtained by generating a synthetic 
companion model and directly extracting the normalized visibilities via Fourier techniques, 
described in the preceding paragraphs. 

We present the azimuthally averaged visibilities from our disk plus companion model 
with best-fit values for companion stellocentric radius and PA in Figure [HI for comparison 
with our traditional disk models (Figures El HI and [6]). We present a synthetic image of this 
model in Figure [9j 

We also examined the closure phases of our data. Closure phases — a measure of the 
complex phase that removes the phase ambiguity produced by the turbulent atmosphere — are 
used to discover and constrain asymmetric structure in observed objects (e.g., a companion). 
We found that the closure phases of our data are insensitive to structures on the small 
scales we are examining. The closure phases are sufficiently noisy to be largely insensitive 
to asymmetric structure more compact than ~100 mas (our putative companion lies at 
~65 mas). 



4. Discussion 



The presence of a compan ion in the TW Hya system is consistent with theories of 



trans ition disk dissipation (e.g., iGoldreich fc Tremaindll982t iBryden et al.lll999t iRice et al. 
20031 ). The fitted stellocentric radius of our putative companion is ~3.5 AU, inside the 
optically thick disk wall at ~3.9 AU. The Keplerian period for such a circular orbit would 
be ~2860 days, assuming a mass for TW Hya of 0.7 M Q (IWebb et al.lll999l ). The time 
between our two epochs is ~720 days, and so during this time the change of the PA of the 
companion would be ~91 degrees. Since we do not have an estimate of PA for our second 
epoch, we cannot test for this orbital evolution in our data. Though the 2009 data do not 
constrain a value for PA, the data are consistent with orbital motion of a companion. This 
is because the companion's flux contribution at 11.66 /zm and 18.30 /zm is less significant; 
another epoch of 8.74 /im observations would enable this exercise. 

We have not quoted uncert ainties on any of ou r fitted parameters. Our parameteriza- 
tions for the reproduction of the iRatzka et al.l (120071 ) model and our new disk model include 
7 and 9 freely varying parameters, respectively. With parameter spaces of increasing dimen- 
sionality, model uniqueness and degeneracy of parameters becomes increasingly problematic. 
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Fig. 8. — This figure shows the synthetic visibilities generated by our disk plus companion 
model, compared to our data. While we fitted models to data in azimuthal wedges, we plot 
the azimuthally averaged data for simplicity. The points with error bars in the figure show the 
normalized, calibrated visibility amplitude as a function of baseline, as in Figures [31 HI and|6j 
As in the previous figures, we show our four separate observations: 8.74 (blue, top left), 
11.66a (black, top right), 11.66b (black, bottom right), and 18.30 (red, bottom left). Our data 
are indicated by unconnected points, and the synthetic model visibilities are shown with 
solid lines. Note that we have not accounted for potential disk variability between the two 
observed epochs, and hence our model tends to fit the average of the two 11.66 /im datasets. 



^ ■ ' I I , ... I ■ ... I I I , ... I ■ ... I I I ■ A 



-5 -4 -3 -2 -1 1 2 3 4 5 
x [AU] 



Fig. 9. — A synthesized image of our disk + companion model at 8.74 /mi. For model 
parameters, see Table O Most of the model components can be clearly seen; only the 
completely cleared region inside the dust sublimation radius (~0.02 AU) is absent at the 
displayed resolution. Several model components are visible: the central star; the optically 
thin disk extending from the dust sublimation radius to the transition disk wall at 3.9 AU; the 
directly illuminated, optically thick wall at 3.9 AU (faint, but visible, at 8.74 /mi), represented 
as a single temperature ring; and the point source companion, producing optically thick 
hot emission at 3.5 AU, required to explain our very resolved 8.74 /im observations. The 
companion has the same surface the hot ring model described in Section 13.21 and 

Table [31 This panel shows our model at 8.74 /mi; the best-fit position angle of the companion 
is shown for this epoch. In order to make the geometrically small elements of the model visible 
in this illustration, we have convolved the model with a PSF with full width half maximum 
(FWHM) ~10 mas. 
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Table 3. Model Properties 









Calvet / Uchida 


Ratzka 


New Disk Model Disk - 


- Companion Model 


(1) 


Kurucz Scaling 




2.0 


1.9 


1.0 


1.0 


(2) 


RL 


(AU) 


0.02 




0.02 


0.02 


(3) 


in 


(K) 


1364 




1364 


1364 


(4) 


S in / S * n , Calvet 




1 




0.79 


0.79 


(5) 


flwall 


(AU) 


3.3 


0.6 


3.9 


3.9 


(6) 


/wall 




0.11 


0.051 


0.17 


0.11 


(7) 




(K) 


130 


271 


116 


116 


(8) 


^atm 


(AU) 




0.6 






(9) 


/ atm 






0.11 






(10) 


T* 

atm 


(K) 




186 






(11) 


Satm / £ m ,Calvct 






1.2 






(12) 


^hot 


(AU) 






2.9 


3.5 


(13) 


/hot 








le-4 




(14) 


Thot 


(K) 






570 


570 


(15) 


Opacities 




C 


R 


C 


C 


(16) 


Position Angle (2007) 


C) 








90 


(17) 


Position Angle (2009) 


(°) 











Note. — The models used to generate SEDs and visibilities presented in this work. Columns separate models that are 
"Calvet-like" , "Ratzka-like" , or newly presented in this work. Rows describe the different components in the model, and 
are (a * indicates the component is optically thin, otherwise it is optically thick): (1) the scaling factor applied to the 
Kurucz model, described in Section[3] (2) the dust sublimation radius, i.e., the inner radius of the optically thin disk; (3) 
the temperature at the dust sub limation radius; (4) the surface density of the inner, optically thin disk, compared to the 
surface density used bv lCalvet e t al. (20q2); (5) the radius of the directly illuminated optically thick wall; (6) the ratio 
of width to radius for the transition disk wall; (7) temperature of the optically thick wall; (8) radius of the optically thin 
wall atmosphere; (9) the ratio of width to radius for optically thin wall atmosphere; (10) temperature of the optically 
thin wal l atmosphere ; (11) the surface density of the optically thin wall atmosphere, compared to the surface density 
used by ICalvet et all (120021 ') ; (12) radius of the optically thick, hotter-than-equilibrium ring inside the transition disk 
wall. This component is only used for the model presented in this work; (13) the ratio of width to radius for the optically 
thick, hotter-than-equilibrium ring; (14) temperature of the optically thick, hotter-than-cquilibrium ring; (15) opacities 
used in each model — "C" for Calvet, "R" for Ratzka. See Scction l3,l.ll fl6) The best-fit position angle of our companion, 
for our 2007 data (17) The best-fit position angle of our companion, for our 2009 data 
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High dimensional chi-squared surfaces are complex, and may contain many local minima. To 
obtain our best-fit parameters, we employ grid-based chi-squared minimizations. We obtain 
a best-fit PA and stellocentric radius of the companion by simply substituting the hot ring in 
our new disk model (§ 13. 2p for the emission from an unresolved source. In doing this, we fix all 
of the many free parameters in our model (radius of optically thick disk rim, temperature of 
model components, etc.). Obtaining the uncertainties for PA and radius from a chi-squared 
surface assuming only one or two varying parameters instead of the actual, larger number, 
would not provide an accurate estimate of the error in these parameters. Since modeling the 
multi-dimensional chi-squared surface is computationally prohibitive using our methods, we 
do not report formal confidence intervals on our model parameters. 

While the flux from the outer regions of the disk (at 11.66 fim and longer w avelengths) is 
rather sensitive to changes in i? W aii and T waU (and has been modeled previously, ICalvet et al. 



(120021 . e.g.)), the quality of the fit to the data depend less strongly on T comp and -R C om P , 



due in part to the variable nature of TW Hya at the shorter wavelengths. We find that 
while changes of just ~10s of Kelvins or ~tenths of an AU in the disk cavity wall can 
lead to an unacceptable fit, changes of a similar magnitude in the hot disk or companion 
cause a much smaller effect. In fact, a temperature for the hot companion as low as 430 K 
produces a lower-quality but still acceptable fit to our data: a deficit of flux at wavelengths 
shorter than 6-8 /im can be explained by a larger stellar contribution. Though some model 
parameters are degenerate, we use this range of temperatures (430-500 K) to investigate 
potential observables predicted by planetary models. 

Even if we assume our lowest acceptable value for T comp , the fitted temperature of 
our companion is hotter than the disk equilibrium temperature at its stellocentric radius 
(~150 K). Given our fitted surface area of emission, and the assumption that this object 
emits like a blackbody, the l uminosity of this object is ~3-6xl0 -4 L & . The age of the 
TW Hya system is ~ 10 Myr JWebb et al.lll999h . which puts an upper limit on the age of 
the putative companion. 

These properties of our proposed companion are consistent with models of planetary 
formation. The temperature of our companion (of age < 10 Myr) is consi stent with the 



predic ted temperatures from planetary thermal evolution models presented i nlFortnev et al. 



(2008) for a planet beginning from the core accretion formation models of iHubickyj et al. 



(120051 ) . The luminosity, however, is sufficiently large to require "hot start" models with 
arbitrary initial conditions; in this case, our compa nion's luminos i ty is consistent a mass of 
8-10 Mj (IFortnev et alihoosl ). Indeed, the work of iMarlev et all (J2003) suggests that such 
a large luminosity is only possible with hot-start models, or if the com panion is undergoing 
an accretion shock phase during the core-accretion formation process. IMarlev et al.l (120071 ) 
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state, however, that luminosities for the hot-start planetary models are highly uncertain and 
depend strongly on the initial entropy of the evolution tracks; this makes mass determination 
of very young planets difficult. We also note that the surface flux density of our companion 
as a function of wa velength is consistent with those of a 500-700 K model presented in 
Fortney et al.l (120081 ) . except for effects due to opacity, which we do not include (see their 
Figure 3). 

The luminosity of our planet is feasible if it is undergoing an accretion shock phase, 
though the m agnitude of an objec t undergoing this phase is highly uncertain, and this phase 
is short-lived (IMarley et al.l 120071 ). Another possible explanation for the high magnitude of 
the companion luminosity is that the companion is surrounded by and accreting circum- 
planetary material; an increase in the surface area of the emitting region would produce a 
larger luminosity. 

Additionally, we have detected significant variability at 11.66 /im in our multi-epoch 
observations. The origin of this variability is unclear. The fitted size (using a ring emitting 
at a single wavelength) for the more resolved epoch is 2.56 ± 0.23 AU, compared to the less 
resolved size of 0.71 ± 0.50 AU. At 11.66 /im, emission is dominated by the optically thin 
transition disk cavity, at relatively small stellocentric radii; indeed, much of the emission 
at this wavelength comes from the hot inner edge of the optically thin disk, at the dust 
sublimation radius. The var i ability is thus likely due to changes in the properties of this 
inner region. iMuzerolle et al.l ( 120091 ) report "remarkable mid-IR variability" in the transition 



disk LRLL 31, in the same wave length range as our observed variable sized emission and on 
timescales as short as one week. iFlaherty fc Muzerolld ( 120101 ) propose models to explain this 
sort of mid-IR variability, using non-axisym metric perturba t ions i n the disk (e.g., a warp 
with va riable scale heig h t, or spiral wave). IFlaherty et al.l (120 111 ) expands on the earlier 



work of IMuzerolle et al.l ( 120091 ) with extensive observations of LRLL 31 over many epochs 
in the infrared. They find that the dust destruction radius stays relatively constant, that 
accretion and IR excess vary over timescales of ~weeks, and that changes in scale height 
and/or warping of the inner disk is likely responsible for the observed variation. They rule 
out accretion and stellar winds as a cause for the disk changes, and conclude that the most 
likely explanation for their observations is a companion or a dynamic interface between the 
stellar magnetic field and the disk. In particular, they describe how a companion, orbiting 
at an inclination relative to the disk, can drag dust from the disk midplane in a periodic 
fashion, producing variable IR emission. Given that the inner disk contributes most of the 
flux at 11.66 fim in our model, variability in the disk itself — rather tha n variable flux from a 
companion — is the most likely explanation, as in IFlaherty et al.l ( 120111 ). As we only have one 
epoch of 8.74 /mi data, we cannot determine if the variability we observe persists at shorter 
wavelengths, though TW Hya has also been observed to be variable in the near-IR and optical 
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as well ( lEisner et al.l I2010L and referenc es therein). A multi-epoch, multi- wavelength study 
of TW Hya, as in iFlaherty et al.l ( 1201 ll ) for LRLL 31, would yield a greater understanding 
of the physical mechanisms responsible for this variability. 

Recently, other compa nions have been detected or inferred around transition disk ob- 
jects. iHuelamo et al.l ( 120111 ) presented evidence of a companion around T Chamaeleontis at 
AL = 5.1 mag and at a separation of 6 2±7 mas using sparse aperture masking and adaptive 
optics at the VLT. Eisner et al.l ( 120091 ) present spatially resolved mid-IR observations from 
T-ReCS of SR 21, at a separation < 100 mas. They predict t hat in the K band the c ompan- 
ion has ~ 5% of the flux of the central star, and ~ 25% at L. iKraus fc Ireland! (120111 ) detect 
a companion around LkCa 15 using non-redundant aperture masking interferometry at the 
Keck-II telescope. This object has projected separations of 72, 101, and 88 mas at thr ee 
different epochs, with a contrast in the L band of AL = 4.7 mag. IKraus fc Ireland! (j201lf )'s 
detection of a companion around LkCa 15 shows evidence of asymmetry and variability in 
stellocentric radius and flux (though they state that some of the asymmetrical smearing is 
indicative of poor data quality). They suggest that circumplanetary material is expected 
around the detected companion, as both planet and star are accreting protoplanetary disk 
mass. 

We predict that in the L band the companion around TW Hya will have ~0. 8-3.0% of 
the flux of the central star (AL ~5.3-3.6 mag), and in the M band (4.7 /xm) the companion 
will have ~2-7% of the flux. Our companion would lie at a separation of ~65 mas from 
the central star. The contrast would be even more favorable at longer wavelengths; the 
hot companion's contribution to the SED peaks at around 7 /iin. Due to its proximity then, 
TW Hya offers the possibility to detect a companion at a smaller stellocentric radius than, but 
at comparable angular separation to, these recent companion discovery publications. High 
resolution mid-IR observations (perhaps speckle interferometry or non-redundant aperture 
interferometry at the Large Binocular Telescope) should verify our predictions. 



We note that lEvans et al.l (120121 ) have placed a lower limit on the //-band contrast of 
a companion in a separation range of 40-80 mas around TW Hya at 5.28 mag, consistent 
with but on the faint end of our L-band contrast estimate of ~5.3-3.6 mag. There are 
several reasons why our prediction for companion contrast at L band could be too large, 
however. For example, if some of the emission from the companion is optically thin, the 
companion flux at L band would be smaller (while maintaining the same flux at 8 jttm), 
leading to a larger value for the contrast at L while preserving the quality of our fit to the 
SED. Some optically thin emission would be evidence that this emission is partially due to 
circumplanetary matter, which would be more spatially extended to produce the same flux 
as a blackbody component. 
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Conclusions 



We present new mid-IR, spatially resolved measurements of the transition disk object 
TW Hya, taken in a novel observing mode at the Gemini telescope using the T-ReCS instru- 
ment. We observed this object using speckle imaging and reduced the data using Fourier 
image analysis techniques. Our individual exposures were short enough to freeze the at- 
mosphere's turbulence, allowing for sub seeing-limited observations. In fact, due to high 
precision calibration of the PSF using our Fourier methods, we probe spatial scales at the 
diffraction limit of the Gemini telescope. At all our observational wavelengths, we resolve 
the science target, TW Hya. 

We recreate and present simple models of TW Hya's disk from the literature. We ana- 
lyze the compatibility of these models with our new data, and show that existing models do 
not reproduce our very resolved 8.74 /im emission. We create a new model that satisfactorily 
explains all the available data: our new resolved mid-I R measurements, long baseline inter- 
ferometric measurements at 7 mm ( IHughes et al.l 120071), the flux density in the mid-IR, and 
long baseline mid-IR visibility amplitudes presented by IRatzka et al.l (120071). This model has 
a lar ge, relatively empty cavity as shown in previous works ( jCalvet et al. 



2002: Uchida et al. 



2004 : IHughes et al.ll2007l ). but also includes a self luminous companion at stellocentric radius 
of ~3.5 AU to explain our hi ghly resolved 8 . 74 Atm data. We note also that the behavior of 
our data is similar to that of iRatzka et al.l (|2007j)'s (8.74 /im emission more resolved than 
11.66 /im emission), who obtain mid-IR observations at another epoch (2005) and at longer 
baselines (~45 m). 

In summary, we present new observations with 8.74 /im emission more resolved than 
11.66 /i m emission at baselin es of ~6 m. This unexpected result is consistent with other 
findings: IRatzka et al.l ( 120071 ) show ~8 /im emission more resolved than ~12 /im emission at 
much longer baselines. If a companion is responsible for this emission, observations should 
show an asymmetry, most significant at the peak emission wavelengths of the companion. 
Our data are consistent with this expectation of asymmetry, but we cannot constrain the 
details of this asymmetry well. We provide estimates for separation and flux ratio of a 
putative companion, but these are uncertain due to large number of components in our 
model, and because the companion is not well-resolved. Evidence for a luminous source at 
a large stellocentric radius is fairly strong, but our constraints on its properties are weaker. 

We measure the size of the emission in our observations at each wavelength, and detect 
temporal variability at the >2.5-<r level between two epochs at 11.66 /im. We speculate that 
this variability reflects variable accretion through the inner disk, signatures of dynamical 
perturbations, and / or perhaps a variable brightness of a self-luminous companion. 
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